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Abstract 

Nuclear interactions of antiprotons in atomic states are discussed. The 
total as well as partial widths for single nucleon capture events are calcu¬ 
lated. These are compared to the X-ray and recent single nucleon capture 
data. The rates of the neutron or proton captures test nuclear density dis¬ 
tributions at the extreme nuclear surface. Recently found cases of neutron 
and proton haloes are analysed. 


1 Introduction 

It has been known for a long time that hadronic atoms are a way to test the nu¬ 
clear surface region: the tail of nuclear density distribution, its isospin structure 
and correlations [1]. There are two methods to learn these properties: 

1) Measurements of the X ray cascade in hadronic atoms that provide atomic 
levels and widths. Some fractions of the level energies are due to nuclear inter¬ 
actions and some parts of the widths are due to nuclear captures of the hadron. 
For highly excited atomic states the nuclear effects are small, for low levels the 
nuclear capture probability increases rapidly with decreasing orbital radii and 
the cascade terminates suddenly. These natural limitations allow to measure 
one level width per atom. Only in some special cases two widths and one level 
shift may be obtained. The levels in question are of large angular momenta 
that locate the nuclear interactions on the nuclear surface. At first, the atomic 
data are used to learn the strength and form of hadron optical potentials. Next, 
some of the level widths may be used to test the nuclear density tail. The shifts 
are usually difficult to interpret and provide a check on the optical potential. 

2) Measurements of the nuclear capture products. A unique detrmination 
of the emitted particles may discriminate captures on protons from captures 
on neutrons and signal nuclear correlations. Many experiments using various 
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detection techniques have presented data, in principle, more informative than 
the X-ray data. Unfortunatly, these are also more difficult to interpret as the 
initial atomic states are not known and the final state interactions, in particular 
the charge exchange reactions, are uncertain. 

Some related results in the p and other atoms are discussed in reviews [2,9] 
and other talks at this meeting. This paper discusses the neutron density dis¬ 
tributions tested by recent antiprotonic CERN experiments, [3]. The latter are 
of the second kind and follow previous antiprotonic [4] and kaonic [5] stud¬ 
ies. Now, however, the separation of pp and pn annihilation modes is done 
in a different way. Instead of the final state mesons it is the final nucleus 
that is detected by radiochemical methods. In this way the rates of reactions 
p{N,Z) mesons{N — 1,Z) and p{N,Z) mesons{N, Z — 1) are found. 
The ratio of these rates allow to detect the number of neutrons relative to the 
number of protons in the region of nuclei where the p capture occurs. The 
measurements done in nuclei where the method is applicable yield results that 
differ widely. For some nuclei (^®Ni,®® Ru) one finds the n/p ratio about unity, 
in heavy nuclei Th,^^® U) it is as large as 5 or 8 while in ^^"‘Sm it 

happens to be significantly less than unity. Once the final and initial states are 
known and the reaction mechanism is understood one can determine where the 
capture occurs and what nuclear region is tested by these experiments. Then 
one can interpret the n/p ratio in terms of ’’neutron halo” or ’’neutron skin ” 
and give more precise meaning to these terms. The purpose of this work is the 
presentation of the basic elements of such an analysis [6]. 

The known difficulties: uncertainty of the capture state and the necessity to 
describe the final state interactions are still present. In particular, the radio¬ 
chemistry is selective, only cold final nuclei are seen and a proper description 
of the final states becomes a question. Fortunately an additional constraint fol¬ 
lows from the measurement itself. It is given by the ratio of two capture rates: 
the rate for single nucleon captures that end with cold residual nuclei and the 
total capture rate. The former, cold captures, ammount to 10-20 percent of 
the total, and are almost independent on the nucleus. If the capture occurs 
from a definite atomic state the total rate is given by the level width and is in 
principle measurable by the X-rays. Such clean experiments are not feasible, 
yet. The chances of nuclear capture from various initial p states have to be 
calculated with an optical potential derived from other atomic data. An exper¬ 
imental check is expected with forthcoming CERN experiments which hopefuly 
will supply also the transition probabilities per single stopped antiproton [6]. 

Next sections describe briefly: the optical potential, final state mesonic in¬ 
teractions and the interpretation of the neutron haloes. 
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2 Nuclear Interactions of Atomic Antiprotons 

Antiprotons captured into atomic orbits cascade down, emitting Auger electrons 
and X-rays, to be finally absorbed by the nucleus. This happens in atomic 
states of high angular momenta, presumably in circular orbits. Such a scenario 
is formed by the atomic cascade that tends to populate states of high angular 
momenta 1. Because of the centrifugal barriers and large NN absorptive cross 
sections, the nuclear interaction of p is rather well localised at distances as 
large as twice the nuclear radius. An important aspect of such peripherality is 
that it allows low density approximations in the theoretical description: quasi- 
free interactions and a single particle picture of the nucleus. It also facilitates 
the description of final mesons, an important issue in understanding of the 
absorption experiments. On the other hand, the surface studies are complicated 
by the sensitivity of results to an uncertain range of hadron-nucleon forces. 

Here, we present a phenomenological description of the antiproton absorp¬ 
tion by nuclei. The level widths are calculated in terms of p optical potentials. 
The simplest one, linear in nuclear density, is of the form [7] 

—* ^TT —♦ 

V°P^{R) = - t^ripiR) ( 1 ) 

^^NN 

where pjsfN is the reduced mass, p{R) is a nuclear density and is a complex 
scattering length. The density p{R) in Eq.(|^) is not the ’’bare” nucleon density 
Po{R) but a folded one 


p{R) = J dupo{R — u)v{u) (2) 

where u is a formfactor that describes the NN force range. The length in 
Eq.® is extracted from the most precise X-ray measurements done in oxygen 
isotopes, [7]. The best fit to these data yields of about —1.5 —i2.5 fm, [7,8]. 
The corresponding potential is deep and black. For the central densities, 
is 200 MeV strong and the mean free path is well below 1 fm. However, both 
the form and the strength of are tested only in the surface region. Thus, 
is determined by the atomic level widths, given by 

r = [ dRp{R) I |2 (3) 

where 'I'jy(i?) is the atomic wave function. Since Tjy R} is determined 
essentially by the angular momentum I and only high momentum states are 
available the absorption strength is peaked at the surface. 

A typical nuclear interaction region in ®®Ni is shown in Fig.l, where the 
absorption density W = p \ '^ \’^ R? \s plotted. There are two special atomic 
states in the capture process. One is called an ’’upper” level which usually is 
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the last one that can be detected by the X rays before the cascading down p is 
absorbed. Its width is obtained from the intensity loss of X-ray transitions. In 
this and in many other nuclei the absorption is most likely to happen from this 
level. The next circular state below it, ’’the lower state”, may be reached in 
some nuclei. In such a case one measures the shape of X-ray line and extracts 
the level width and shift. Nuclear absorptions may happen also in higher atomic 
orbits in a way that is not detected by X-ray studies. The chances for the p to 
reach the low levels in question are not known well. 



R (fm) 


Figure 1: The antiproton absorption densities from n=6, 1=5 (upper level) in 
®®Ni: Ws for the NN annihilation range Vo = .75 fm, p is a ’’bare” neutron 
density. The cold antiproton absorption densities on a neutron (integrand in 
Eq.(§): Ai for the NN annihilation range ro = 1 fm and Ag for the range 
To = .75 fm. Normalisations are arbitrary. Missing probabilities (left scale): 
Tmiss continuous is due to phase space alone, Pmiss dash-dotted is calculated 
with corrections for the experimental pion momentum distribution. The flat 
dashed curve is Pdh from the HF model. 

The peripherality of capture depends on the range of NN forces. The range 
parameters in Eq. i) may be adjusted to fit the atomic and low energy scatter¬ 
ing data. Gaussian profile formfactors exp(—(r/co)^) have been used, [8], and 
typical best fit values are: rd ss 1 fm (for SE) and Cor « 1.5 fm (for 5RE). On 
the other hand, calculations based on the NN potentials yield values Coi of 0.7 
fm up to 1.5 fm, [11], for different NN states and different ways to go off-shell. 
The latter values are probably the upper and lower limits of Coi, while the best 
fit number is located in-between. An effect of the range uncertainty is shown in 
Fig.l for a partial decay width. The effect on the full width is essentially the 
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same. 

Optical model calculations based on the NN interaction potentials [10,11] 
show that the lengths are not the NN S-wave scattering lengths. The 
latter are smaller and repulsive (positive ). The effective is density de¬ 

pendent and has a complicated structure. At the nuclear surface it reflects a 
long attractive tail of the pion exchange forces, at distances around the nuclear 
radius it may turn to repulsion due to repulsive scattering lengths, and it is 
rather uncertain at the nuclear matter densities. On the other hand, the phe¬ 
nomenological best fit represents cumulative effect of the S and P wave 

absorptive amplitudes and can be well understood in terms of the free At tv at • 
Theoretical optical potentials are more complicated than formula (|^, but do 
not reproduce the data as accurately as the latter with the best fit parameters. 

The level widths reflect all modes of nuclear absorption. The initial stage of 
this process, an elementary NN annihilation, generates w 2 GeV energy that 
turns mostly into kinetic energy of the final state mesons. These mesons excite 
residual nuclei by inelastic processes. To calculate the widths one sums over 
unobserved nuclear excited states. The large energy release allows a closure 
over these states and the effective is close to the absorptive part of the 

free This is not true when final nuclear states are limited by the mea¬ 

surements. Thus, the radiochemical detection allows only cold nuclei i.e. nuclei 
either in the ground states or excited less than the neutron separation threshold, 
[3]. The effect of this limitation is discussed in the next section. 

3 Nuclear NN Annihilation and Final State In¬ 
teractions 

The aim of this section is to calculate the rate of nuclear p annihilations that 
lead to cold final nuclei. This is done in few steps: 

1) An amplitude for the NN annihilation into mesons is assumed 

and introduced into the nuclear transition amplitude in the impulse approxima¬ 
tion. 

2) The emission probabilities are calculated and summed over mesonic and 

nuclear final states. For an isolated NN annihilation this procedure would 
produce the absorptive cross section and via unitarity condition the absorptive 
amplitude AtTVTv- I^or nuclear captures leading to cold nuclei we limit the sum¬ 
mation over final states to the states of elastic meson nucleus scattering. This 
limited summation generates the again, but now it is folded over nuclear 

final state interaction factors, Pmiss(^), which describe the probability that the 
annihilation mesons born at point R miss the residual nucleus. 

3) Let a nucleon N occupy a single particle level a with a wave function 
(pa{X). Then, the pN annihilation process limited by Pmiss leads to a single hole 
final nuclear state. The experimental condition allows only those initial levels 
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a that end in cold final nuclei. A factor Pdh(-A) = ‘Pa accounts 

for this condition. 

4) The finite range effects are described by separable potentials as done in 

[ 11 ]- 

The result for a partial width corresponding to a cold capture on a nucleon 
s = (proton or neutron) is now 

r,(co/d) = [ I |2 V{Y - X)ps{X)PMX)PmUR) (4) 

Pnn J 

where is an effective length and R is the birthplace of the mesons. The 
assumption, justified later, is that all the mesons are emitted from the central 
point of the annihilation region R = . 

Eq.(|^ is the result, now we explain briefly the final state interactions that 
determine Pmiss and in the next section we turn to calculations of the nuclear 
densities. 

The spectrum of mesons consists essentially of pions correlated in a sizable 
fraction into p and uj. These heavy mesons propagate some 1 fm and then turn 
into pions. The pion multiplicities range from 2 to 8 with an average 4—5 and 
an average momentum as large as 400 MeV. Nuclear interactions of these pions 
may be absorptive, inelastic or elastic. All the absorptive and almost all the 
inelastic processes would not leave the residual nuclei in cold states and so the 
production rate for the latter is given essentially by the elastic scattering. This 
allows an optical potential description. In addition, in the bulk of phase space 
the pions are fast enough to allow also an eikonal description. Following this, 
the wave function for each pion is taken in the form 

= exp(ipif- iS{p, 0) (5) 

where p is a momentum, ^ is a coordinate and S is expressed in terms of the 
pion-nucleus optical potential 

pOO 

S{p,0 = / ds{ 

Jo 

The S is calculated by integrating the local momentum over the stright line 
trajectory. Due to nuclear excitations and pion absorptions this wave is damped 
with a rate described by 3S'. The latter is generated by absorptive part of the 
pionic optical potential This damping follows the whole path but the 

main effect comes from regions of large nuclear densities and not the region 
around the birth place We assume that all functions S{p,^) are related to 
the central point of annihilation R that is the NN CM coordinate. With the 
mesonic wave functions d) the dependence on the total momentum of mesons P 
factorizes approximately to a plane wave form. One consequence is that the NN 
CM ’’conservation” S{R—R') arises in the transition probabilities integrated over 


(p^ — C/°P* i^+ ps) — p) 


( 6 ) 
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final momenta. Now, the final state pion interaction factors may be collected 
into a probability distribution 

Pmiss(-R) =< n I exp(-5(K, R)) |2> (7) 

It is a product of the eikonal factors averaged over the number and phase space 
of final pions with some allowance for an unknown momentum dependence gen¬ 
erated by 

The calculations of Pmiss are performed in a Monte Carlo procedure. The 
optical potential for pions must cover wide momentum range from the threshold 
up to 0.9 GeV but the phase space favours a region just above A resonance. This 
potential is related to the pion nucleon forward scattering amplitudes and in this 
way to the pion nucleon cross sections. That method is established around the 
A [13]. Here, it is extended to higher N^i, resonances which are described 
by Breit-Wigner amplitudes. The two nucleon absorption mode is taken in a 
phenomenological form [14]. Performing these calculations one finds that: high 
energy expansion of the square root in Eq. is satisfactory, higher resonances 
cannot be neglected and the black sphere limit is a good approximation in dense 
regions. The result is close to a pure geometrical estimate that relates Pmiss(P) 
to the solid angle of the nucleus viewed from the point R [9]. In the surface 
region of interest the Pmiss is a linear function of the radius. It is very fortunate, 
it makes this calculation fairly independent on the size of the NN annihilation 
region and on the uncertain range of the heavy meson propagation. 

Examples of Pmiss, Pdh and the density A for cold absorption generated by 
Eq. are given in Fig.l. The latter is seen to be more peripheral than the total 
absorption density W. The uncertainty due to the NN force range is rather 
small, beeing moderated by effects of the strong absorption. The uncertainty 
of the initial atomic state is shown in Table 1. The ratio of cold to total cap¬ 
ture rates T'^/T* raises quickly with the increasing angular momentum. The 
experimental results, given in Table 2, indicate that a significant contribution 
of high I states is unlikely. On the other hand, some participation of lower I 
states is possible. However, due to the atomic cascade properties, that seems to 
be unlikely. A definitive answer requires experimental and theoretical studies 
which are beeing undertaken, [6]. 

The absorbtion widths are given by superpositions of high moments of the 
nuclear density distributions. We find the 21 — 2 moment as the dominant one 
for the total widths and the 21 moment to dominate the cold capture width. 
The ’’neutron halo ” is thus understood as a ratio of these high moments of 
neutron and proton density distributions. The ’’neutron skin ” is related the 
m.sq. radius or other low moments. 
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Table 1: Column 2 contains main quantum number n for the lower and upper 
circular states. The weigths give probabilities for nuclear capture calculated 
under the assumption that the circular atomic state n = Uupper + 1 is fully 
occupied. The absorption widths: total T* and cold T'^ are calculated with the 
AD model, Rnp = .63. 


ELEMENT 

n 

weight 

r7r‘ 

Tn/Fp 

58Ni 

4 

0 

.095 

.69 


5 

.16 

.097 

.69 


6 

.83 

.11 

.70 


7 

.01 

.15 

.71 


8 

0 

.22 

.71 

90 zr 

6 

.24 

.106 

4.67 


7 

.72 

.128 

5.30 

i^^Sm 

7 

.01 

.087 

3.65 


8 

.75 

.099 

3.98 

238]j 

9 

.29 

.106 

6.55 


10 

.71 

.138 

8.24 


4 Nuclear Densities 

As the simplest estimate we use an asymptotic density (AD) model. It follows 
the Bethe-Siemens approach [12] although the larger input includes: charge 
density distribution, neutron and proton separation energies and difference of 
the rms radii of proton and neutron density distributions. At central densities 
a Fermi gas of protons and neutrons is assumed. The Fermi momenta are 
determined by the densities and Fermi energies are fixed by the separation 
energies. This gives depth of the potential well which in the surface region is 
extrapolated down in the Woods-Saxon form. The densities are given by the 
exponential damping of the nucleon wave functions due to the potential barriers. 
For protons a Coulomb barrier is added and potential parameters (c, t) are fitted 
to reproduce the experimental charge density down to 5 percent of the central 
density. For neutrons the same t is used but c is chosen to obtain the rms 
radius equal (or larger by 0.05 fm in the heaviest nuclei) to the proton density 
rms radius. This model is expected to generate average level densities, it misses 
shell effects and correlations. 

A second model used to determine neutron and proton densities is the 
Hartree-Fock (HF) and the Hartree-Fock-Bogolyubov (HFB) scheme with the 
effective two-body Skyrme-type interaction. Our aim in using HF and HFB 
methods to find nucleon densities at the extreme tails of the nuclear matter 




distribution (at distances of 8—15 fm from the center) is rather unusual. The 
necessary practical condition is the use of a HF code not restricting the asymp¬ 
totic form of s.p. wave functions. This excludes e.g. all codes using the harmonic 
oscillator basis. In the present work we have applied the code solving HF equa¬ 
tions on the spatial mesh, in which all fields and densities are expressed in the 
coordinate representation. 

The HF method disregards completely pairing correlations. To account for 
them we used the HFB theory [16,17] which unifies the self-consistent description 
of nuclear orbitals, characteristic of HF method, and the mean field treatment of 
residual pairing interaction into a single variational theory. The effective force is 
the ten-parameter Skyrme SkP interaction described in ref [18]. It has a virtue 
that the pairing matrix elements are determined by the force itself, contrary to 
other Skyrme-type interactions which define only the particle-hole channel. The 
paired HFB ground-state has not the BCS form so there is no simple pairing 
gap parameter although a kind of average gap can be defined as the pairing 
potential average over occupied states. 

The most severe restriction of the presented results is the imposed spherical 
symmetry, both in the HF and HFB codes. It allows enormous simplification of 
solutions, in particular the HFB equation takes the form of two coupled differ¬ 
ential equations in the radial variable for each value of s.p. angular momentum. 

The density matrix is obtained by summing contributions from the lowest 
s.p. orbits. The degeneracy of the spherical subshells is handled by taking 
contribution of the last orbit in the filling approximation, i.e., an appropriate 
occupation probability smaller than one is, if necessary, associated with this 
orbit. 

Some results are collected in Table 2. To find F^/Fp a ratio of pn and pp ab¬ 
sorptive amplitudes Rnp is needed. One number Rnp = .63 follows from mesonic 
studies in [4]. It depends on uncertain final charge exchange processes and 
we use the deuteron value Rnp = .82 [15]. The latter coincides with a good fit 
to F„/Fp in 5®Ni. 

Our conclusions are: 

1) The AD model overestimates the neutron haloes. The latter are not 
given by the binding energies and Coulomb barriers alone. The shell effects 
(angular momentum) are essential, correlations of HFB type have rather small 
effect. These conclusions are similar to the results of subcoulomb neutron pickup 
studies [19]. 

2) Even at these nuclear peripheries at least 2-3 nucleon orbitals are involved 
in the capture. 

3) There are two anomalous cases: Yb and Te. The anomalies are apparently 
related to a strong E2 mixing in those atoms. 

4) An interesting case of proton halo is seen in It is not understood 

yet. 
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Table 2: Experimental and calculated results for and r„/rp. Calculations 
are averaged over few atomic orbitals weighted as in Table 1, with Rnp = .82. 



pc/pt 

Pn/Pp 

pc/pt 

Pn/Pp 

pc/pt 

Pn/Pp 


Exp. [3] 


Asympt. 


HP 


58Ni 

.098(8) 

0.9(1) 

.11 

.90 

.110 

.785 

90 zr 

.161(22) 

2.6(3) 

.12 

4.9 

.125 

2.54 

96Ru 

.113(17) 

0.8(3) 

.10 

1.7 

.099 

.944 

i30Te 

.184(36) 

4.1(1) 

.12 

2.6 

.124 

3.14 

144 Sm 

.117(20) 

< .4 

.09 

1.9 

.094 

1.38 

i54Sm 

.121(20) 

2.0(3) 

.10 

5.1 

.110 

3.34 

176Yb 

.241(40) 

8.1(7) 

.12 

4.8 

.111 

3.23 

232Th 

.095(14) 

5.4(8) 

.12 

7.6 

.087 

3.80 

238-u 

.114(9) 

6.0(8) 

.13 

10 

.092 

4.09 


5) The antiprotonic study of nuclear surface is a promising method, despite 
some uncertainties. 
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